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The problem of noise-induced escape from a metastable state arises in physics, 
chemistry, biology, systems engineering, and other areas. The problem is well un- 
derstood when the underlying dynamics of the system obey detailed balance. When 
this assumption fails many of the results of classical transition-rate theory no longer 
apply, and no general method exists for computing the weak-noise asymptotics of 
fundamental quantities such as the mean escape time. In this paper we present a 
general technique for analysing the weak-noise limit of a wide range of stochastically 
perturbed continuous-time nonlinear dynamical systems. We simplify the original 
problem, which involves solving a partial differential equation, into one in which only 
ordinary differential equations need be solved. This allows us to resolve some old 
issues for the case when detailed balance holds. When it does not hold, we show 
how the formula for the mean escape time asymptotics depends on the dynamics of 
the system along the most probable escape path. We also present new results on 
short-time behavior and discuss the possibility of focusing along the escape path. 
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I. INTRODUCTION 



The phenomenon of escape from a locally stable equilibrium state arises in a multitude of 
scientific contexts If a nonlinear system is subjected to continual random perturbations 

('noise'), eventually a sufficiently large fluctuation will drive it over an intervening barrier 
to a new equilibrium state. The mean amount of time required for this to occur typically 
grows exponentially as the strength of the random perturbations tends to zero. 

Research on this phenomenon has focused on the case when the nonlinear dynamics of 
the system in the absence of random perturbations are specified by a potential function. 
In a recent paper 0] we have introduced a new technique for computing the weak-noise 
asymptotics of the mean first passage time (MFPT) to the barrier. Our technique, unlike 
the bulk of earlier work, is not restricted to the case when the zero-noise dynamics arise 
from a potential. Because of this we can readily and quantitatively treat systems without 
'detailed balance,' whose dynamics are determined by non-gradient drift fields, or are oth- 
erwise time- irreversible. We deal here with overdamped systems, in which inertia plays no 
role. Overdamped systems without detailed balance arise in the theory of glasses and other 
disordered materials , chemical reactions far from equilibrium || , stochastically modelled 
computer networks [0-|9|], evolutionary biology [|K|, and theoretical ecology fllTf . 

In the multidimensional escape problems most frequently considered in the literature, 
the most probable escape path (MPEP) in the limit of weak noise passes over a hyperbolic 
equilibrium point (i.e., saddle point) of the unperturbed dynamics. In the case of non- 
gradient drift fields exit through an unstable equilibrium point can also occur [|J. Other 
new possibilities, such as exit through a limit cycle HT2|Jl3|] , arise as well. However, if the 
unperturbed dynamics are determined by a potential, exit in the limit of weak noise must 
occur over a hyperbolic equilibrium point, and the asymptotics of the MFPT are given by 
a classic formula, originally derived in the context of chemical reactions by Eyring |14 . 



Over the years the Eyring formula has been rederived and generalized by a variety of 
alternative approaches p5HT7|. To illustrate its use, consider a two-dimensional system 
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whose dynamics are specified by a sufficiently smooth drift field u = u(x, y) symmetric 
about the X~£IX1S clS displayed in Fig. [I]. S — (x$, 0) and H = (0, 0) denote the stable and 
hyperbolic equilibrium points, and the barrier lies along the |/-axis. The position of a point 
particle representing the system state, moving in this drift field and subjected to additive 
white noise w(t), satisfies the Ito stochastic differential equation |TB| 



dxi(t) = Ui(x(t)) dt + e 1//2 <Tj dwi(t), % — x,y. 



(1) 



Here a x and a y quantify the response of the particle to the perturbations in the x and y- 
directions; the corresponding diffusion tensor D is diag(_D x , D y ) = diagfa^ 2 , a y 2 ) and will in 
general be anisotropic. If the drift field u is obtained from a potential function = <p(x, y) 
by the formula Ui = ~Didi<p then u will in general not be a gradient field, but the relation 



between u and D will ensure detailed balance 19 



If detailed balance holds, the Eyring formula for the e — > asymptotics of the MFPT r 



is 
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where (du y /dy)(S) is du y /dy evaluated at the stable point S, and (du y /dy)(H) is du y /dy 
evaluated at the hyperbolic point H, etc. In this example the MPEP lies along the x-axis; 
this fact was used implicitly in the expression for the Arrhenius (i.e., exponential) factor 
in (0). If the vector field were asymmetric and the MPEP were curved, the Arrhenius factor 
would involve the integral J(D _1 u) • <ix taken along the MPEP [^0 . 

Eq. (||) displays a frequently encountered feature of weak-noise escape problems: the 
dependence on u of the pre-exponential factor in the MFPT asymptotics is limited to a 
dependence on the derivatives of u near the stable and hyperbolic equilibrium points. If de- 
tailed balance is absent this will generally not be the case pl| , p2| . However, most quantita- 
tive work on the case when detailed balance is absent has dealt with one- dimensional state 



spaces, due to the difficulty of higher- dimensional computations. Talkner and Hanggi |23 



computed MFPT asymptotics for a multidimensional model of an optically bistable system 
without detailed balance, in which the drift field did not arise from a potential. But due 
to the simplicity of their model, the prefactor did not display a complicated dependence on 
the drift. Multidimensional models without detailed balance have not been treated in full 
generality. 

In Sections |1| and [ill] we start from scratch and compute weak-noise MFPT asymptotics 



for multidimensional systems by singular perturbation methods, using matched asymptotic 



expansions. Similar approaches have been used before |12|J24|1 , but our treatment yields a 
systematic method of computing the pre-exponential factor, even in cases when detailed bal- 
ance is absent. As noted, we have previously applied our techniques to the nonclassical case 
of exit over an unstable equilibrium point (i.e., an equilibrium point at which the lineariza- 
tion of u has only positive eigenvalues) ||]. The MFPT asymptotics differed considerably 
from the classical escape formulae such as Eq. (|2|); the prefactor in the asymptotics of r de- 
pended on e. In this paper we flesh out the ideas in our previous paper by applying them to 
drift fields with the standard structure of Figure [TJ. We discover that even when the MPEP 
passes over a saddle point, if detailed balance is absent the MFPT asymptotics likewise 
differ from those given by the Eyring formula. Although the prefactor is independent of e, 
it depends on the behavior of u along the entire MPEP. We quantify this dependence; it can 
be much more complicated than a dependence on the derivatives of u at the endpoints. 

Our treatment reveals, however, that when u is symmetric about the MPEP as displayed 
in Fig. m the prefactor can readily be computed by integrating two coupled ordinary differen- 
tial equations along the MPEP. Previous work has left the impression that partial differential 
equations must be solved. In the absence of detailed balance, the MFPT asymptotics are 
determined by a quasipotential |25[ , the solution of a nonlinear partial differential equation. 
However the derivatives of this function along the MPEP, which turn out to determine the 
prefactor, satisfy ordinary differential equations. 

The prefactor is very sensitive to the behavior of the drift field in a neighborhood of the 
MPEP. Integrating the ordinary differential equations along the MPEP requires knowledge of 



the extent of 'differential shearing' along the MPEP, i.e., knowledge of the second derivatives 
of u there. In many cases differential shearing can give rise to a focusing singularity along 
the MPEP, in which case more sophisticated techniques must be used. 

When detailed balance is absent, in the limit of weak noise the short-term dynamics 
of the system, as well as the MFPT asymptotics, display novel behavior. As noted, the 
system fluctuates around the locally stable equilibrium state many times before it finally 
escapes along the MPEP. If the drift field u arises from a potential, each of the unsuccess- 
ful escape attempts follows an outward 'instanton' trajectory anti-parallel to some integral 
curve of u and then falls back along the same integral curve (an 'anti-instanton' trajectory). 
In the absence of detailed balance, the time-irreversibility will generally cause the instanton 
trajectories not to be anti-parallel to the anti-instanton trajectories, which remain integral 
curves of the drift field. Hence unsuccessful escape attempts will proceed with high proba- 
bility along closed loops containing nonzero area, in contrast to the more familiar situation 
for systems with detailed balance. We discuss this in Section [F7J conclusions appear in 
Section [V]. 



We now derive the weak-noise MFPT asymptotics for any u with the general structure of 
Fig. |l], but not necessarily derived from a potential. We allow anisotropic diffusion, so long 
as the principal axes of the diffusion tensor are aligned with the MPEP. After approximating 
the probability density of the system along the MPEP, we will compute r by the Kramers 
method of computing the probability flux through the separatrix, i.e., the boundary of the 
basin of attraction of the stable point. 

The Fokker-Planck equation for the probability density p is p = C*p, with 



II. THE ANALYSIS 



C* = (e/2)(D x d x 2 + D y d y 2 )-u x (x,y)d : 



•x 



u y (x,y)d l 
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(d x u x (x,y)) - (d y Uy(x,y)). 



(3) 
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If absorbing boundary conditions are imposed on the separatrix (in this case, the y-axis), the 
MFPT in the weak-noise limit may be approximated by (— Ai) -1 , where Ai is the eigenvalue 
of the slowest decaying density mode p\ (i.e., the eigenfunction of C* whose eigenvalue has 
the greatest real part). Ai is negative and converges to zero as e — > 0, so in the weak-noise 
limit (— Ai) -1 displays Arrhenius growth. We have 

/oo / roo poo 

(D x e/2)d xPl (0,y)dy/ / Pl (x,y)dydx (4) 
-oo / JO J—oo 

since the right-hand side is the normalized flux of probability through the separatrix. The 



leading asymptotics of Ai, as computed by this formula, will be unaffected |12] if p\ is taken 
to satisfy £*pi = 0- 

According to (Q), to approximate Ai we must approximate the normal derivative of p\ 
along the separatrix. But since the MPEP passes through a hyperbolic point, the probabil- 
ity density will be concentrated in a small region (of size 0(e 1//2 )) about this point as e — > 0. 
It therefore suffices to approximate p\ near the hyperbolic point. To evaluate the denomi- 
nator of (^), we must also approximate p\ near the stable point. The latter approximation 
is straightforward: near (xs,0) we take 

Pifay) ( 5 ) 
~ expHAT^A^Kx - x s ) 2 + D y - l \\ y {S)\y 2 )/e] 

where we now write X X (S), X y (S), X X (H), and X y (H) for the (real) eigenvalues of the drift 
field linearized at the stable and hyperbolic points, respectively. This Gaussian approxima- 



tion permits the evaluation of the denominator of (|). It is 7rey D x D y / y X x (S)X y (S) + o(e), 
since contributions from other regions are negligible. 

It remains to approximate p\ in the vicinity of the hyperbolic point. As a necessary first 
step, we approximate p\ along the MPEP, i.e., the x-axis. In the vicinity of the MPEP we 
use a standard WKB approximation 

pi(x,y) ~ K(x,y)exp (-W(x,y)je) . (6) 
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Here W is the quasipotential, the 'nonequilibrium potential' investigated by Graham and 
others |26| , |27j| . K satisfies a transport equation, and W an eikonal (Hamilton- Jacobi) equa- 
tion: H(x, VW) = 0, with H the Wentzell-Freidlin Hamiltonian [E5f 



tf(x, P ) = ^ 2 + ^ 2 + u(x). p . (7) 

So W(x,y) can be viewed as the classical action of the zero-energy trajectory from (xs, 0) 
to (x,y). If detailed balance holds, i.e., Ui = —Didi<j) for some potential function 0, it is 
easily checked that W(x, y) = 2<ft(x, y) + C, for C an arbitrary constant. If detailed balance 
is absent then W will in general be more difficult to compute. The zero-energy classical 
trajectories determined by the Hamiltonian of (^) include the instanton trajectories; we will 



return to this point in Section IV 



Given the drift field structure shown in Fig. |], and assuming sufficient smoothness, we can 
expand u near the x-axis in powers of y: 

u x = v (x)+v 2 (x)y 2 + O(y 4 ) (8) 
u y = Ul (x)y + 0(y 3 ). (9) 

Using this notation, we have, e.g., X y (S) = Ui(xs) and X X (H) = dv /dx(0). The function t> 2 
measures the extent of differential shearing along the MPEP. 

Again assuming sufficient smoothness, we can expand W itself near the x-axis in powers 
of y: 

W(x, y) = f (x) + f 2 (x)y 2 + 0(y 4 ) (10) 

where f 2 {x) is a measure of the transverse behavior of the 'WKB tube' of probability current 
along the MPEP; clearly, /^(^O = \d 2 W / dy 2 (x , 0). Physically, f 2 measures the transverse 
lengthscale on which probability density is nonnegligible; the tube profile will be approxi- 
mately Gaussian, with variance proportional to e/f 2 (x) at position x. 

Substituting the Ansatz (g) into C*p± = and equating the coefficients of powers of e 
yields a system of equations for the the functions f , f 2 , K in terms of v , v 2 , U\. (The system 



closes: no higher derivatives of u or W enter.) The first of these is the one-dimensional 
eikonal equation H(x, /o(x)) = 0, where 

H(x,p x ) = ^y?x 2 + v (x)p x (11) 

is a Hamiltonian governing motion along the MPEP. For this Hamiltonian, 

x = D x p x + v (x) (12) 

is the expression relating velocity and momentum. 

By the eikonal equation, fo(x) can be viewed as the momentum p x of an on- axis classical 
trajectory with zero energy. It follows from the Hamiltonian flllD that there are only two 
solutions to the eikonal equation: Jq(x) = —2vq{x)/D x (arising from an instanton trajectory 
moving against the drift, with x = —Vq(x)), and / = (arising from an anti-instanton 
trajectory following the drift, with x = +Vq(x)). That these trajectories are directly anti- 
parallel to each other is a consequence of our assumption that the drift field is symmetric 
about the x-axis MPEP. Since physical exit trajectories begin at the stable point and ter- 
minate near the hyperbolic point, when constructing the WKB approximation (||) we use 
the instanton solution /q(x) = —2vq{x)/D x rather than the anti-instanton solution /o = 0. 

The equation for f 2 is a nonlinear Riccati equation, and may be written as 

f 2 = -2D y f 2 2 - 2uJ 2 + 2v v 2 /D x (13) 

where the dot signifies derivative with respect to instanton transit time. (Since the instanton 
trajectory satisfies x = —v (x), the left-hand side equals —Vof 2 .) This equation is new, and 
has an immediate physical interpretation: it describes how the WKB tube spreads out or 
contracts under the influence of its environment, as one progresses along the MPEP. The 
final inhomogeneous term on the right-hand side of (|T3| ) is particularly important. If the 
differential shearing v 2 is sufficiently negative in some portion of the MPEP, this term can 
drive f 2 to zero before the hyperbolic point is reached 0,0. If this occurs, to leading order 
the WKB tube splays out to infinite width. The tube width actually remains finite, but 
becomes larger than 0(e 1 ^ 2 ); this becomes clear if higher-order terms are taken into account. 
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If goes negative, further integration of (|T3|) will normally drive f% to — oo in finite 
time; for a pictorial example of this 'focusing' phenomenon, which is really the appearance 
of a singularity in the nonequilibrium potential, see Fig. 1 of Day |28[]. Eq. ( [3~3|) accordingly 
gives a new quantitative measure of the validity of our WKB approximation: for it to be 
valid, f2 must remain positive along the entire MPEP. 

It has not been recognized that the validity of the WKB approximation along a straight 
MPEP can be so easily checked. Nonlinear equations resembling ([i~3| ) have been derived by 



Schuss [|18] and Talkner and Ryter |13|j29|l , but in the context of motion along the separatrix 



rather than along the MPEP. Their equations are homogeneous rather than inhomogeneous, 



so they give rise to no singularities. Ludwig and Mangel |2Jj3^] derived a homogeneous 
equation in the context of motion along an anti-instanton trajectory, i.e., motion following 
the drift. 

The third and final equation is for the function K. The transport equation for K(x, y) 



is a well-known partial differential equation fl2|,[T3||, but for the computation of the MFPT 



asymptotics an ordinary differential equation along the MPEP will suffice. We may approx- 
imate K within the WKB tube as a function of x alone, in which case we find 

K/K = - Ul - D y f 2 (14) 

In the vicinity of the MPEP, the WKB approximation @ arising from perturbations of the 
on-axis instanton trajectory x = —Vq(x) is completely determined by the ordinary differen- 
tial equations ( pT3[) and ([14]), together with the initial conditions /^(^s) = \\ y (S)\/D y and 
K(xs) = 1- These follow from the requirement that the WKB approximation match the 
Gaussian approximation ([|) near the stable point (x$,0)- 

Finally we can approximate p\ near the hyperbolic point H = (0,0). In the diffusion- 
dominated region within an 0(e l l 2 ) distance of the origin, pi satisfies the equation 

(e/2)(D x d 2 xPl + D y d 2 yPl ) - d x (\ x x Pl ) - d y {\ y y Pl ) = (15) 

which is separable: P i(x,y) = P *(x) P \(y). By the symmetry of u, P \ must be even; the 
absorbing boundary condition on the separatrix implies that P \ must be odd. The separation 
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constant, and the overall normalization, are found by requiring that this solution match the 
WKB approximation (|6]). 

For the matching to occur, the separation constant must equal zero. The symmetry 
requirements mandate that 



Pl (x,y) ~ Ce x2 ' A y2{\,X) exp {-\\{H)\y 2 / D v 



(16) 

We have introduced the rescaled variable X = xJ2\ x (H)/D x e, and 2/2(^5 ") is the odd 
parabolic cylinder function |31] of index |. 2/2(^5 *) can be expressed in terms of elementary 
functions, but we write p\ in terms of it to facilitate comparison with our earlier work on exit 
over an unstable equilibrium point B. There a similar expression occurred, but the index of 



the parabolic cylinder function was not fixed. Cf . the treatment of Caroli et al. Jl5| , in whose 
treatment Weber functions (i.e., parabolic cylinder functions) of variable index were used. 

72 e x2 / 4 as X -> +00. 

+00 asymptotics of 



Ttl 



Note that y' 2 (^,0) = 1 by definition, and that y 2 (^,X) 

The normalization constant C is obtained by matching the X 
px(x,y) to the WKB solution near the hyperbolic point. This gives 

2 



C = K(0)J2/tt exp 



u (x) dx 



(17) 



D x e Jo 

where K(0) must be computed by integrating ([14]) along the MPEP from S to H. Since f 2 
appears on the right-hand side of (|T4]) , this in turn requires an integration of (13) along the 
MPEP. 

The essential point here is that the behavior of p\ near the hyperbolic point — in par- 
ticular, its overall normalization — can in general only be obtained by integrating along the 
MPEP, and will depend on the entire history of u and its first and second partial derivatives 
along it. Substituting our approximations (H) and (|TB|) for p\ into Eq. (^) now yields 



-~—J\K(S)\\ x (H)< 

T 7T v 



\ \{H) 



x exp 



D; 



e jo 



uq(x) dx 



(18) 



which resembles the Eyring formula (0) but contains a new 'frequency factor' K(0). 
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Formulae similar to Eq. fll8|) appear in Schuss |18| and Talkner |L3]. However, Eq. (|18| 



does not indicate how Po(H) is to be found. Our treatment makes it clear that for drift 
fields with the structure of Fig. [I], the frequency factor K(0) and the extent to which the 
MFPT asymptotics differ from the Eyring formula are most easily computed by integrating 
the ordinary differential equations flT3| ) and ( fn| ) along the MPEP. 

III. EXPLICIT EXAMPLES 

If the drift field is derived from a potential, we recover the Eyring formula as follows. 
As noted, if Ui = —Didi4> f° r some potential function then W(x, y) = 2<f)(x,y) up to an 
additive constant. This simplifies the calculation of f2~- the nonlinear Riccati equation for f 2 
need not be integrated explicitly, though it could be. Necessarily 

d 2 <p{x,y) 



, ld 2 W(x,y) 



2 dy 2 
-D 



X dUy 



v 



dy 



(x,0) 

u^/Dy. (19) 



(*,0) 



Eq. QTSJ) settles a recurring issue. It is sometimes assumed [HJ , in studying multidimensional 
escapes, that if u varies too rapidly along the MPEP the WKB approximation may break 
down due to f 2 going negative at some point along the MPEP. We have just shown that if 
detailed balance holds, this will never happen, so long as U\ satisfies the obvious transverse 
stability condition of being strictly negative. In the language of chemical physics, that 
f 2 (x) = —U\[x)jD y for every x says that in this case, transverse fluctuations are in local 
thermal equilibrium at all points x along the MPEP. If detailed balance is lacking, fi{x) 
will not depend on U\{x) alone; Eq. (|13|) makes this observation quantitative. 

It follows immediately from ( |T3D and ([19]) that K = constant; since K{x$) = 1, A"(0) = 1 
also. So the formula (|18|) reduces to the Eyring formula when u is derived from a potential. 

To illustrate the power of our technique, we now examine an overdamped system without 
detailed balance where the frequency factor K(Q) cannot be computed analytically. Suppose 
for simplicity that D x = D y = 1, and consider the drift field 
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u x = x — x 3 — axy 2 (20) 
u y = -y- x 2 y. (21) 

It is easily checked that this has the structure shown in Fig. p] for any a, with xs = 1, and 
that for a = 1 (and only a = 1) this u is derivable from a potential <fi. In this case 

0(x,y) = -ix 2 + ix 4 + i 2 / 2 + |xV. (22) 

Note that for all a, the eigenvalues X X (S) = —2, X y (S) = —2, X X (H) = 1, and X y (H) = —1. 

If a = 0, is independent of y and the escape time problem becomes essentially 
one-dimensional. In particular, the MFPT asymptotics are given by the one-dimensional 
Kramers formula [|J 

' ' '\X X (S)\X X (H) exp(--A^). (23) 



r 7r " \ e 



As a varies between zero and one, we expect the asymptotics of r to interpolate smoothly 
between those of the Kramers formula (|23|) and those of the Eyring formula (|2]). 

Our technique can be used to solve for K(0) and the prefactor in the MFPT asymptotics 
for any value of a. In general, a numerical integration of the equations (|T3| ) and ([14]) for f'2 



and K is required. Though this is straightforward, it is more illustrative of the value of our 
technique to perturb about the values of a at which analytic solutions can be obtained. 

We therefore consider the drift field (|20|),(|2lD with a — 1 + 8, |5| <C 1, and solve for the 
first-order (in 5) correction to K(0) in ([18]). We have u± = —(1 + x 2 ), v = x — x 3 , and 
V2 = — (1 + S)x. Assuming sufficient smoothness, we expand /2 in powers of 5: 

h = /i 0) + 5/ 2 (1) + 0(6) 

= ( 1 + x i) + SfP + (6) (24) 

since f% = —U\/D y when a = 1. Substituting this expansion into (0), we find a linear 
equation for the first-order correction f% 1 . 

^/« = 2(i±4)/« + 2x (25) 

\ A, 1 Ay' 1 A.' 1 A.' 
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This inhomogeneous first-order equation is easily solved: 



= TT^W [l°g* + (1 - * 2 ) + i(z 4 - 1)1 • (26) 



ar 



Notice that when 5^0 and detailed balance is lacking, to first order f2{x) will differ 
from —u\(x)/D y at all x along the MPEP. So if q ^ 1, we expect that the transverse 
fluctuations are in local thermal equilibrium nowhere along the MPEP. 

Using Eq. (11), and converting the derivative with respect to instanton transit time to a 
space derivative gives the equation for K: 

K'/K = 5— (log x + § + K - x 2 ) . (27) 

Its solution is 

K(x)/K(l) (28) 



exp 



-S(l - x 2 - (2a; 4 - 4a; 2 ) logx)/4(l - x 2 f 



From this we find, to order S, K(0) = K(l)(l + U) = 1 + §5. So 



1 1 + 15 



T IT 



y/\K(S)MH) 



Xy(S) expf--A0) (29) 



€ 



to first order in 5. Eq. (^9| ) displays the first-order correction to the Eyring asymptotics. 

Just as we have perturbed about the analytically soluble a = 1 model, we could also 
perturb about the analytically soluble a = model, where detailed balance is altogether 
absent. When a = the nonequilibrium potential W(x, y) cannot be computed explicitly. 
However its transverse second derivative fi along the MPEP can be computed explicitly by 
integrating ([13]). The computation of the first-order correction f% yields corrections to the 
Kramers asymptotics of (|23|); details are left to the reader. 

Notice that since t>2 = —<yx, when a is sufficiently large and positive the function f 2 
obtained by integrating the differential equation ( |T3D along the MPEP will develop a zero 
at some point between x — 1 and x = 0. At this point the WKB tube will (naively, see 
Section splay out to infinite width, and the WKB approximation will break down. This 
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is qualitatively reasonable: by (0), if a > the differential shearing near the MPEP is 
such as to enhance the likelihood of motion away from the MPEP in order to overcome the 
resistance of the drift field. 

IV. SHORT TIME BEHAVIOR AND UNSUCCESSFUL ESCAPES 

We turn to the short-term dynamics of our model, and consider the trajectories followed 
by the particle in the exponentially many unsuccessful escape attempts which precede the 
final successful escape along the MPEP. 

The situation in the presence of detailed balance is well known. An unsuccessful escape 
attempt, defined as a sample path which leaves the diffusion-dominated region of size 0(e 1 ^ 2 ) 
surrounding the stable point S, with high probability follows a trajectory moving anti- 
parallel to the drift. The return path follows a deterministic trajectory satisfying x = u(x). 
This symmetry is a consequence of time-reversibility. 

If detailed balance is absent, the instanton trajectories (i.e., the classical trajectories with 
nonzero momentum and zero energy, as determined by the Wentzell-Freidlin Hamiltonian, 
emanating from the stable point) serve as the most likely unsuccessful escape paths. This 
can be seen in several ways. The Wentzell-Freidlin Hamiltonian fl7|) is dual to the Onsager- 
Machlup Lagrangian 

= ~ Ux \ 2 + ^D~\y ~ u y\ 2 ■ ( 30 ) 

The probability that a sample path will be close to any specified trajectory x*(t), < t < T, 
should to leading order fall off exponentially in e, with rate constant equal to the integral of 
the Lagrangian along x*. This integral is minimized when x* is a classical trajectory, and 
it is easily seen that minimizing the integral over transit time T, for fixed endpoints x*(0) 
and x*(T), selects out the classical trajectories of zero energy. 

There are in general many such trajectories, including the deterministic trajectories which 
simply follow the drift. (For them, L is identically zero.) We have adopted the field-theoretic 
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nomenclature of Coleman [^], according to which the most likely exiting trajectories are 



called 'instantons,' and the return trajectories are anti-instantons. A careful treatment, using 
functional integral methods, shows that the most likely trajectories, i.e., the local minima 
of the action functional, consist of any number of (pieces of) instanton trajectories, each 
followed by the corresponding anti-instanton trajectory: the relaxation toward the stable 
point S provided by the deterministic dynamics. Our nonequilibrium potential W(x,y) is 
simply the classical action of the instanton trajectory terminating at (x,y). 

Graham and others [2f| have stressed the importance of the nonequilibrium potential in 



determining the long-term dynamics or MFPT, as well as the MPEP. However its impor- 
tance in short-term dynamics has seldom if ever been clearly enunciated. If the system is 
observed to be in some internal state S' far from the stable state S , with high probability it 
reached S' by following an instanton trajectory, and will return to S by following an anti- 
instanton trajectory. This applies to states S' which are outside the diffusion-dominated 
zone of size 0(e 1 ^ 2 ) surrounding S, and the phrase "with high probability" can be given a 
quantitative meaning in the weak-noise limit. In the absence of detailed balance the outward 
instanton trajectories are not anti-parallel to the inward anti-instanton trajectories, and the 
resulting closed loops contain nonzero area. This is a consequence of time-irreversibility. 

Fig. |2| illustrates this phenomenon. If D x = D y = 1 and the drift field u is given by ( |20|) 
and fl2lD with a = 0, detailed balance is absent. Though the nonequilibrium potential cannot 
be computed explicitly except along the x-axis, the instanton and anti-instanton trajectories 
are easily computed numerically. Typical trajectories are shown. That the x-axis instanton 
trajectory (i.e., the MPEP) is atypical, being directly anti-parallel to the drift, is in part a 
consequence of the symmetry of our model. It is also a consequence of the fact that by (§) 
and (0), we are assuming du y /dx = du x /dy to hold along the MPEP. Within the WKB tube 
of width 0(e 1/2 ) about the MPEP, to high accuracy u is irrotational and detailed balance 
holds. However the small deviations from detailed balance due to differential shearing can, 
as we have seen, prevent transverse fluctuations from being in local thermal equilibrium. 

Fig. ^| also illustrates the phenomenon of focusing: several of the instanton trajectories 
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displayed intersect each other, though not along the MPEP. Due to such intersections (which 
typically form caustics p3| , |34]j , or singular surfaces of codimension unity) the computation 
of the nonequilibrium potential W will in general require a minimization over trajectories 
which are piecewise classical (with zero energy) rather than classical. Some preliminary 
investigations of this effect have been made by Graham and Tel [3"5||, but it is not yet 



clear how to handle foci appearing along the MPEP. In our model it is easily checked that 
/2 being driven through zero to — oo is the sign of a focus along the MPEP, in fact of a cusp 
catastrophe ||33|| . The appropriate extensions to the WKB approximation are now under 



study. 

Note that though our instanton and anti-instanton trajectories are incident on the stable 
point S and/or the hyperbolic point H, the portion of any trajectory within an 0{e 1 ^ 2 ) 
distance of either equilibrium point is without physical significance. On this lengthscale 
diffusion dominates in the sense that the particle can diffuse from any point to any other 
within 0(1) time. 

This sheds light on the seeming paradox of the trajectories having infinite transit time. 
Since x = u(x) for anti-instanton trajectories an infinite amount of time is needed to reach S\ 
similarly an infinite amount of time is needed for the instanton trajectories to emerge from 5". 
The physical escape trajectories, both successful and unsuccessful, should be viewed as 
emerging from the diffusion-dominated zone rather than from S itself. So for example any 
unsuccessful escape trajectory directed along the x-axis will have transit time 

|A x (S)|- 1 log(l/ e 1 / 2 )+0(l). (31) 

The transit time of the return path, i.e., the corresponding anti-instanton trajectory, will 
be the same. Similarly the amount of time needed for the final, successful escape trajectory 
to approach H will be 

A x (i/)- 1 log(l/ e 1 / 2 ) + 0(l). (32) 

In all 
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P^r 1 + \ X {H)- 1 ] logil/e 1 / 2 ) + 0(1) (33) 

time units will be required for the MPEP to be traversed in full during the successful escape. 

We see that there are two timescales: the exponentially large MFPT and the logarithmic 
timescale on which escape attempts occur. In the weak-noise limit the latter is very brief 
compared to the former; the number of unsuccessful escape attempts grows exponentially. 
The successful escape, when it finally occurs, will on the MFPT timescale occur almost 
instantaneously; this justifies the term 'instanton.' 



V. SUMMARY AND CONCLUSIONS 

We have seen that even in the absence of detailed balance, the weak-noise MFPT asymp- 
totics of a nonlinear system with a point attractor can readily be computed. There is no 
simple expression for the pre-exponential factor in the asymptotics, such as that given by 
the Eyring formula. Rather, it follows by integrating the differential equations (0) and ( |Hj) 



along the MPEP. This technique clarifies the new phenomena that can occur, such as the 
formation of focusing singularities due to differential shearing along the MPEP. The effects 
of these singularities on the MFPT asymptotics will be treated in a future publication. 

Our treatment makes it clear that the time-irreversibility present in systems without 
detailed balance manifests itself even unto the shortest time scales. It is well known that 
in such systems stationary states, and the quasi-stationary state we employ to compute the 
MFPT asymptotics, are characterized by a global circulation of probability. However the 
circulation discussed in Section |TV| differs from the conventional sort |3(J in that it occurs 



locally, at the level of individual unsuccessful escape attempts. Recognition of this point is 
essential for a proper computation of the MFPT asymptotics. 

In this paper we have treated the two-dimensional case on account of its simplicity, but 
our treatment generalizes to higher dimensions. Curved MPEPs, diffusion tensors whose 
principal axes are not aligned with the MPEP, and nonconstant diffusion tensors can all be 
treated by appropriate extensions of the techniques presented here. In higher-dimensional 
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models, or models with curved MPEPs, the Riccati equation (p~3|) must be replaced by a 
matrix Riccati equation. This extension is related to a standard result of stochastic control 
theory |J7|: the effects of perturbing about optimal trajectories, obtained by solving a 
Hamilton- Jacobi equation, can be computed by solving a matrix Riccati equation. 

Most studies of the overdamped escape problem have dealt with time-reversible stochastic 
models, with detailed balance. This is a reflection of the origins of the escape problem 
in statistical mechanics, where nonlinear dynamics are usually derived from a potential. 
In broader applications of stochastic modelling, which include engineering as well as the 
biological and social sciences, there is little reason to expect time-reversibility. Our approach 
should prove useful in the wider context. 
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FIGURES 

FIG. 1. A drift field u symmetric about the x-axis, with stable equilibrium point S = (xg,0) 
and hyperbolic point H = (0,0). \ X (S) and \ y (S), the eigenvalues of the linearization of u at S, 
are taken to be real and negative. This sketch assumes X X (S) = X y (S), but that is not assumed in 
the text. 

FIG. 2. Typical instanton and anti-instanton trajectories for the drift field specified by equa- 



tions (|20| ) and (21), with a = 0. The two families are superimposed: the anti-instanton trajectories 
are largely straight, but the instanton trajectories curve strongly to the right. The focus mentioned 
in the text occurs to the right. Diffusion here is isotropic; D x = D y = 1. 
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